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Abstract 

The magnetic moments of heavy Sg baryons containing a single 
charm or bottom quark are calculated in the framework of light cone 
QCD sum rules method. A comparison of our results with the predic- 
tions of other approches, such as relativistic and nonrelativistic quark 
models, hyper central model, Chiral perturbation theory, soliton and 
skyrmion models is presented. 



PACS: ll.55.Hx, 13.40.Em, 14.20.Lq, 14.20.Mr 



*e-niail: taliev@metu.edu.tr 
'I'e-niail:el46342@nietu. edu.tr 
■I- e-niail:ozpineci@nietu. edu.tr 



1 Introduction 



During the last few years, exciting experimental results are obtained in the 
baryon sector containing a single b-quark. The CDF Collaboration observed 
the states and E*^ [Ij, while both DO [2] and CDF [3] Collaborations 
have seen Sf,. Recently, BaBar Collaboration reported the discovery of fi* 
with mass splitting m^, - m^, = (70.8 ± 1.0 ± l.l)MeV 

The masses of the heavy baryons have been studied in the framework 
of various phenomenological models [5]- [13] and also in the framework of 
QCD sum rules method [14j- [26j. Along with their masses, another static 
parameter of the heavy baryons is their magnetic moment. Study of the mag- 
netic moments can give valuable information about the internal structures 
of hadrons. 

The magnetic moments of heavy baryons have been studied in the frame- 
work of different methods. In [271 [25] the magnetic moments of charmed 
baryons are calculated within naive quark model. In [291 EQ], magnetic mo- 
ments of charmed and bottom baryons are analyzed in quark model and in 
[HT] heavy baryon magnetic moments are studied in bound state approach. 
Magnetic moments of heavy baryons are calculated in the relativistic three- 
quark model [32], hyper central model [33], Chiral perturbation model [M] . 
soliton model [25] , skyrmion model [SB] and nonrelativistic constituent quark 
model with light and strange qq pairs [37] . In [38] the magnetic moments 
of Sc and Ac baryons are calculated in QCD sum rules in external electro- 
magnetic field. In [SESni, the light cone QCD sum rules method is apphed 
to study the magnetic moments of the Aq, {Q = c, b) and SqAq transitions 
(more about this method can be found in [HI |12], |l3l HI] and references 
therein ). 

The aim of the present work is the calculation of the magnetic moments 
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of the Sf, baryons recently observed by DO and CDF Collaborations within 
the light cone QCD sum rules framework. The plan of the paper is as follows. 
In section 2, using the general form of the the baryon current, the light cone 
QCD sum rules for and baryons are calculated. In section 3 we present 
our numerical calculations on the and Sc baryons. In this section we also 
present a comparison of our results with the predictions of other approaches. 

2 Light cone QCD sum rules for the mag- 
netic moments 

In order to calculate the magnetic moments of Hg {Q = b, c) in the framework 
of the light cone QCD sum rules, we need the expression for the interpolating 
current of Eg. To construct it, we follow [TT], i.e. we assume that the strange 
and light quarks (sq) in Eq are in a relative spin zero state (scalar or pseudo 
scalar diquarks). Therefore, the most general current without derivatives and 
with the quantum numbers of Sq can be constructed from the combination 
of aforementioned scalar or pseudoscalar diquarks in the following way 

VQ = eatc[{q''^Cs'h, + /?(g«^C75s'')]Q^ (1) 

here a, b and c are color indices, C is the charge conjugation operator, Q = b, 
or c, g = u, or d and (3 is an arbitrary parameter. Having the explicit ex- 
pression for the interpolating current, our next task is to construct light cone 
QCD sum rules for the magnetic moments of Hq baryons. It is constructed 
from the following correlation function: 

n(p, q)=tj d'xe'P^ < 7 I T{r^Q{x)r]Q{0) |}0 > . (2) 

The calculation of the phenomenological side at the hadronic level pro- 
ceeds by inserting into the correlation function a complete set of hadronic 
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states with the quantum numbers of Sg. We get 



< -QAP2j I -qAPiJ >7 



2 2 
P2 - "^L 



2 2 



(3) 



Isolating the ground state's contributions, Eq. ([3]) can be written as 



n 



2 2 
^2 ~ 



< Sq(p2) I Sq(pi) 



< Sq(pi) ho I > 



2 2 



< /li(j92) I KKVx) >7 



(4) 



where pi = p + g, P2 = P and q is the photon momentum. The second term in 
Eq. (jlj) describes the higher resonances and continuum contributions. The 
couphng of the interpolating current with the baryons Sq is determined as 



(5) 



where u^q{p) is a spinor describing the baryon Hg with four momentum p 
and \q is the corresponding residue. 

The last step for obtaining the expression for the physical part of the 
correlator function is to write down the matrix element < Sq(p2) | Sq(Pi) >7 
in terms of the form factors. Using Lorentz covariance, this matrix element 
can be written as 



< ^q{Pi) I 2q(P2) 



^^USgiPl) 



2m- 



f2 



(/l + /2)7m + 



(jh±P2)f^j^ 



U3Q{p2)e'', 

(6) 



where fi{q^) and f2{(f') are the form factors and is the photon polarization 
vector. 
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For calculation of the Sq magnetic moments, the values of the form factors 
only at = are needed because the photon is real in our problem. Using 
Eqs. dUlH]) for physical part of the correlator and summing over the spins of 
initial and final Sg baryons, the correlation function becomes 



n = -Xle^"^ 



2 2 

P2 — 



^1 + 



2m=„ " p? — ml 



(7) 



From this expression, we see that there are various structures which can 
be chosen for studying the magnetic moments of Sq. In the present work 
following pS], we choose the structure ^2^4 that contains magnetic form 
factor /i + /2 and at = it gives the magnetic moment of Hg in units of 
eh/2m'=Q. Choosing this structure in the physical part of the correlator, for 
the magnetic moments of Hg we obtain 



n = -A^ ^ ^ 1 



(8) 



where /ih^ = (/i + \q2=o are the magnetic moments of Sg in units of 
efi/2mHQ. 

In order to calculate the magnetic moments of Hg baryons, the expression 
of the theoretical part of the correlation function is needed. After simple 
calculations for the theoretical part of the correlation function in QCD we 
obtain 

+ /?2^^^'Tr(75^f 75^r')} I 0) (9) 
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where S'^ = CSjC and C and T are the charge conjugation and transposition 
operators, respectively and Sq and Sq(^s) are the heavy and hght (strange) 
quark propagators. 

The correlation function from QCD part receives three different contribu- 
tions: a) perturbative contributions, b) nonperturbative contributions, where 
photon is emitted from the freely propagating quark (in other words at short 
distance) c) nonperturbative contributions, when photon is radiated at long 
distances. To obtain the expression for the contribution from the emission 
of photon at short distances, the following procedure can be used: Each one 
of the quarks can emit the photon, and hence each term in Eq. ^ corre- 
sponds to three terms in which the propagator of the photon emitting quark 
is replaced by: 

St d'yF^'^y.Sf^^^ix - y)i,Sf^^^{y)^ ^ , (10) 

where the Fock-Schwinger gauge, x^A^{x) = has been used. Note that the 
explicit expressions of free light and heavy quark propagators in x represen- 
tation are: 

ree _ ^ '^Q 



St 



2 



c/ree rn^ Ki^rnqV^^ ^ 



where Ki are Bessel functions, rriu^d = and rris ^ 0. The expression for the 
nonperturbative contributions to the correlation function can be obtained 
from Eq. ([9]) by replacing one of the light quark propagator by 

s% - -^rr,g^(r,)., , (12) 
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where Tj = |l, 75, 70, 27570, aa/s/V^^ and sum over Tj is implied, and 
the other two propagators are the full propagators involving both perturba- 
tive and nonperturbative contributions. In order to calculate the correlation 
function from QCD side, we need the explicit expressions of the heavy and 
light quark propagators in the presence of external field. 

The light cone expansion of the propagator in external field is obtained 
in [l3]. It receives contributions from various qGq, qGGq, qqqq nonlocal 
operators , where G is the gluon field strength tensor. In this work, we 
consider operators with only one gluon field and contributions coming from 
three particle nonlocal operators and neglect terms with two gluons qGGq, 
and four quarks qqqq [B]. In this approximation the expressions for the 
heavy and light quark propagators are: 



( 2 ,2N2^^''(^^)^/^'^ 



+ 2 ^ .. vx^G^""^^ 
vrir, — k^ 



.1 r ^ 



-igs J du 
m, 



2 oG,,uiux)a^^ - ux^'G^^{ux)Y-r^r^ 



(13) 



where A is the energy cut off separating perturbative and nonperturbative 
domains. 

In order to calculate the theoretical part, from Eqs. (!9l)- (fT3!) it follows 
that the matrix elements of nonlocal operators qTiq between the photon 
and vacuum states are needed, i.e. (7(g) | q{xi)Tiq{x2) \ 0). These matrix 
elements can be expanded near the light cone = in terms of the photon 
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distribution amplitudes 



(7(g)|g(x)a/,i.g(0)|0) = -iegqq{e^q^ - e^q^) J due''"''' yXVy{u) + -j^Mu] 

\ 1^ J Jo 

7(g)|g(x)(7,G^,(t;x)g(0)|0) = -ie,{qq) (e^q, - e,q^) j Pa,e^(°'-+^"«)^^5(a,) 
7(g)|g(x)(7,G^,275(^x)g(0)|0) = -ie,{qq) {e^q, - e,q^) J I)a,e'("'?+''"«)'^^5(«,) 
liqMx)g,G^,ivxh^^,qiO)\0) = e,h^q^ie^q,-e,q^) J Pa,e*("'-+^"«)«M(a,) 
l{qMx)gsG,,{vx)t^^q{0)\0) = e,h,q^{e^q, - e,q^) j P«,e^("'-+''"«)'?^V(«,) 
'^{q)\q{x)ao,i3gsG^uiyx)q{^)\Q) = eg{qq) | (^e^, - (^9au - ^(lo'Xu + Qu^a) 



+ 



qx J 



H iqaXf3 - qpXa){,e^q 

qx 



^uqf,] 



where x is the magnetic susceptibihty of the quarks, ^-yiu) is the leading twist 
2, ip'"{u), ip°'{u), A and V are the twist 3 and h^{u), A, % {i = 1, 2, 3, 4) are 
the twist 4 photon distribution amphtudes (DA's), respectively. The explicit 
expressions of DA's are presented in numerical analysis section. 

The theoretical part of the correlation function can be obtained in terms 
of QCD parameters by substituting photon DA's and expressions for heavy 
and light quarks propagators in to Eq. Q. Sum rules for the Eq magnetic 
moments are obtained by equating two representations of the correlation 
function. The higher states and continuum contributions are modeled using 
the hadron-quark duality. Applying double Borel transformations on the 
variables pi = {p + q)"^ and = both sides of the correlator, for the 

Sq baryon magnetic moments we get: 



+ 



iP' - l)e 



so 



c^'b p[s)ds 



-I 



288772 



7i?(6eQ + e^rrisml <qq> 



{(3' - l)e "'b m\ 



Q 



72M% 



[cs + eg)mo < ss >< qq > r]i 



+ 



432M| 



- l)ml < ss >< qq > 



36eQ + (e^ + eg)A(wo) 



+ {p^ + l)/3^egm,m2 < ss > {r]2 + ^"{uq)) 
+ i - P^) < ss >< qq > 



72 



Ueg - (e, + eg)[?7i - mQXi^yiuo)] 
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(6eQ + es)'lEml < qq > 



- l)e ^'^M^m, 



"■Q 



e s 77i 



A2 

3(/?^ - l)eQm^ < gg ><( -3 + 2-fE + '^H^] 



2887r2 

+ e^mQ < ss > (1 + + CgmQ < qq > 
deQirngrris 



(l-/?^)(7,, + H^]) 
Q 



+ 



1447r2 



< ss > (1 + + 2 < gg > (1 - 



(15) 



where = and Mq = jgr^^. Since the masses of the initial and 

final baryons are the same, we will set Mf = M| and uq = 1/2. The functions 
appearing in Eq. (fTSll are defined as: 



Pis) 



J Vai J dvS{ai)5{aq + vag — Uq), 



- 1) 



144^2^2 \ ^o(6eQ + e,)m, < gg > ln( — 



+ 



+ 



3(l + /?2)eom^ [13 



647r4 



DO TTlf-) O 



(1 — f3'^)ms < qq > ip 
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8 + /n(- 



A2 



2887r2m^ 



- l)ml{es + eeg) < gg > 



+ 6eQ 



3m2 {(1 + < > +2(1 - Z^^) < gg >}^ 



10 
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+ 



'cs + I2ec 



-J — mgrns < qq> 



(1 + 7b)(^22 + 2^ 



12 j 



^02 - V^32 - V^22 - —^"20 + 3(2^32 + ?>i)22 + 



A2 

/n( — 2" 



+ I2eq{^{(3'^ -l)ms < qq> r]ii)2i 



+ (l+/5')/37^2 



2 H 



1 1 m 

^21 - ^10 + 2^20 + 2^00 + 



where, 



da,. 



'0 ^0 

and functions ipnm are defined as 



dag6{l 



a 



all 



2 \n—m ' 

q) 



(17) 



;i8) 



(16) 



Note that the contributions of the terms ~< > are also calculated, but 
their numerical values are very small and therefore for customary in Eq. (fTS!) 
these terms are omitted. From Eq. f|T5|) it follows that for the determination 
of the Eq baryon magnetic moments, we need to know the residue Xq. The 
residue can be obtained from the two-point sum rules and is calculated in 
|25j . For the current given in eq. ([1]) it takes the following form: 

(1 + /3W 



+ 



+ 



e <3 



'■/M 



so 



ds e-^/^^^ 



2H 



X 



X 



--- + 1 



12\nx 



;i + 6')B 



— — 1 — a; 1 + 5x 
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+ 



2%2 



'S 



(1 + l3'^)ml < ss > 
6 



e— _ (1 _ f3^)ml <qq> e~-^y^''l 



+ 



Jo 




(gg)B 



6 



) 



where, x = rriQ/s. 

3 Numerical analysis 

The present section is devoted to the numerical analysis of the magnetic 
moments of Sg baryons. The values of the input parameters, appearing 
in the sum rules expression for magnetic moments are: {uu){l GeV) = 
{dd){l GeV) = -(0.243)3 GeV^ 6(1 GeV) = 0.8{uu){l GeV), ml{l GeV) = 
0.8 GeV^ HZ], , A = 300 MeV and /s^ = -0.0039 GeV^ [16]. The value 
of the magnetic susceptibility x(l GeV) = —3.15 ± 0.3 GeV"^ was obtained 
by a combination of the local duality approach and QCD sum rules |16]. 
Recently, from the analysis of radiative heavy meson decay, x(l GeV) = 
— (2.85 ± 0.5) GeV~^ was obtained ^48j, which is in good agreement with 
the instanton liquid model prediction [19], but slightly below the QCD sum 
rules prediction [16]. Note that firstly the magnetic susceptibility in the 
framework of QCD sum rules is calculated in [50] and it is obtained that 
x(l GeV) = —4.4 GeV~^. In the numerical analysis, we have used all three 
value of X existing in the litrature and obtained that the values of the mag- 
netic moments of baryons are practically insensitive to the value of x- 
The photon DA's entering the sum rules for the magnetic moments of Sg 
are calculated in [16] and their expressions are: 
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ij^iu) = 3 (3(2n - 1)2 - 1) + ^ {I5w}^ - 5w:f ) (3 - 30(2m - 1)^ + 35(2m 



(1 - (2n - Vf) (5(2n - 1)^ - l) ^ (l + ^< - 



^(a-) = 360agaga^ ^1 + w^^{7ag - 3)^ , 

V(aj) = 540w^(ag — aq)aqaqa'^g, 

h^{u) = -10{1 + 2k+)cI{u-u), 

A{u) = AOu'^u^ {3k - k+ + l) 

+8(C2+ - 3C2) [uu{2 + 13mm) 

+ 2^3(10 - 15m + 6u^) \n{u) + 2u^{10 - 15m + 6u^) ln{u)] , 

Ti{ai) = -120(3(2 + C^)(ag - ag)agagag, 

r2(a,) = 30a2(a.-a,)((fi:-K+) + (Ci-O(l-2a,) + C2(3-4a,)), 

%{ai) = -120(3(2 - (2") (ttg - ag)agagag, 

%iai) = 30«J(«5-«,)((fi: + fi:+) + (Ci + Ci^)(l-2ag)+C2(3-4«3)), 

S{ai) = 30al{{K + K+){l-ag) + {Ci + C+){l-ag){l-2ag) 

+ C2[3{aq - aqf - ag{l - ag)]}, 

S{ai) = -30al{{K-K+){l-ag) + {Ci-C+){l-ag){l-2ag) 

+ C2[^{ag- aqf - ag{l- ag)]}. 

The constants appearing in the wave functions are given as [16] </52(l GeV) = 
0, w}^ = 3.8 ± 1.8, = -2.1 ± 1.0, K = 0.2, k+ = 0, Ci = 0.4, C2 = 0.3, 
Ct = and C2+ = 0. 

From the exphcit expressions of the magnetic moments of Hg baryons, it 
follows that it contains three auxiliary parameters: Borel mass squared M^, 
continuum threshold Sq and (3 which enters the expression of the interpolating 
current for Sq. The physical quantity, magnetic moment yUsg, should be 
independent of these auxiliary parameters. In other words we should find 
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the "working regions" of these auxihary parameters, where the magnetic 
moments are independent of them. 

The value of the continuum threshold is fixed from the analysis of the 
two- point sum rules, where the mass and residue Ah^ of the Sg baryons 
are determined [25j, which leads to the value sq = 6.5^ GeV^ for and 
So = 3.0^ GeV^ for Sc. If we choose the value sq = 6.4^ GeV"^ for Eb and 
So = 8 GeV^ for H^. the results remain practically unchanged. Next, we 
try to find the working region of where independent of it at 

fixed value of /3 and the above mentioned values of so- The upper bound 
of is obtained requiring that the continuum contribution should be less 
than the contribution of the first resonance. The lower bound of is 
determined by requiring that the highest power of be less than 30°/o 

of the highest power of M^. These two conditions are both satisfied in the 
region 15 GeV^ < < 20 GeV^ and 5 GeV^ < M| < 8 GeV^ for and 
Sc, respectively. 

In Figs. 1 and 2, we depict the dependence of /i=o and on at fixed 
value of P and So = 6.5^ GeV^. In Figs. 3 and 4, we present the dependence 
of /iso and fi^+ on at fixed value of (3 and so = 3.0^ GeV^. From these 
figures, we see that the values of the magnetic moments of Ef, and E^ exhibit 
good stability when M|, varies in the region 15 GeV^ < < 20 GeV^ 
and 5 GeV^ < M| < 8 GeV^, respectively. The last step of our analysis is 
the determination of the working region for the auxiliary parameter /3. For 
this aim, in Figs. 5, 6, 7, and 8 we present the dependence of the magnetic 
moments of Eq baryons on cos^ where tan^ = (3, using the values of 
from the "working region" which we already determined and at fixed values 

of Sq. 

From these figures we obtained that the prediction of the magnetic mo- 
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ment fi^^ if^Sc) is practically independent of the value of the auxiliary param- 
eter (3. From all these analysis we deduce the final results for the magnetic 
moments in Table 1 for x = —3.15 GeV^. Comparison of our results on the 
magnetic moments of Sg baryons with predictions of other approaches, as 
we already noted, is also presented in Tab lei. 













Our results 


-0.045 ± 0.005 


-0.08 ±0.02 


0.35 ±0.05 


0.50 ±0.05 


RQM [32J 


-0.06 


-0.06 


0.39 


0.41 


NQM [32J 


-0.06 


-0.06 


0.37 


0.37 


[33] 






-1.02^-1.06 


0.45 ^ 0.48 


[34J 






0.32 


0.42 


[35J 






0.38 


0.38 


[36] 






0.28 


0.28 


[37J 






0.28^0.34 


0.39^0.46 



Table 1: Results for the magnetic moments of Sg baryons in different ap- 
proaches. 

We see that within errors our predictions on the magnetic moments are in 
good agreement with the quark model predictions. Our results on the mag- 
netic moments of Sc are also close to the predictions of the other approaches 
except the prediction of [33] on /i=o. 

In summary, the magnetic moments of Hg baryons, which were discovered 
recently (more precisely 3^ was discovered) are calculated in framework of 
light cone QCD sum rules. Our results on magnetic moments are close to 
the predictions of the other approaches existing in the literature. 
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Figure 1: The dependence of magnetic moment /x=;o on at sq = 6.5^ GeV^ 
and P — ±5, — 1. 
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Figure 2: The same as Fig. 1 but for /i^-- 
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Figure 3: The same as Fig. 1 but for /xso and at Sq — 3.0^ GeV^ . 
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Figure 4: The same as Fig. 3 but for 
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Figure 5: The dependence of the magnetic moment /x^o on cos6 at Sq — 
6.5^ GeV^ and for M% = 15 GeV^ and M| = 20 GeV^. 
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Figure 6: The same as Fig. 5 but for 
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Figure 7: The same as Fig. 5 but for /i^o and Sq = 3.0^ GeV'^ and for 
M| = 5 GeV^ and M| = 8 Ge^^. 
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Figure 8: The same as Fig. 7 but for //^-f 
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